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SUMMARY 


This report considers some aspects of the prediction of flap-lag 
instabilities of hingeless rotor blades in hovering flight. Of partic- 
ular interest is the sensitivity of analytical predictions of flap-lag 
stability to various analytical modelling assumptions. The dependence 
of the characteristic modes of motion on blade pitch angle is examined 
for various analytical models using root locus techniques. 

Prediction of flap- lag stability using a single bending mode for 
each degree -of- freedom is examined in the case in which the bending 
modes are assumed to be the same in the flap and lag directions and 
are independent of pitch angle and stiffness distribution. It is 
shown that this model gives results analogous to those obtained by 
Ormiston employing a rigid blade model with the blade and hub stiffness 
represented by springs in the limiting cases of the elastic coupling 
parameter R = 0 and 1. For Intermediate values of R the results are 
shown to be quite different. The mode shape assumptions are shown 
to result in what is referred to as the parallel spring model in 
contrast to Ormiston ’s model which is referred to as a series spring 
model. The similarities and differences between these two models is 
developed in some detail. The differences between these two models 
are examined for various typical rotor blade characteristics. Other 
aspects of the sensitivity of this problem are also considered. 

The notation used and the basic development of the rigid blade 
equations follows that of two papers: "Linear Flap-Lag Dynamics of 
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Hingeless Helicopter Rotor Blades in Hover", by Ormiston and Hodges, 
Journal of the American Helicopter Society, April 1972, and "a Study 
of Stall-Induced Flap- Lag Instability of Hingeless Rotors", by 
Ormiston and Bousman, Journal of the American Helicopter Society, 
January 1975. 

This report also indicates how equivalent lumped spring constants 
and masses are determined from a modal analyses of a rotor blade. 
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INTRODUCTION 


It has been, shown that hingeless rotor blades may experience 
instabilities owing to the coupling between out-of -plane (flap) bending 
and in-plane (lag) bending. The coupling between the flap motion and 
the lag motion arises from aerodynamic, elastic and inertial forces 
acting on the rotor blade. 

The stability boundaries of this coupled motion appear to exhibit 
a considerable sensitivity to the approximations employed in developing 
the equations of motion. This is a' result of the fact that the damping 
of the in-plane bending motion is very small in the uncoupled case and 
therefore, comparatively small coupling terms can destabilize the coupled 
motion. 

This report examines some aspects of the sensitivity of flap-lag 
motion stability to various assumptions which may be employed. 

Of particular interest here is the examination of various modelling 
assumptions associated with describing the structural properties of the 
blade/hub system. 

In order to be able to conduct parametric studies of flap-lag 
stability, one particularly convenient approximation to the bending 
motion of the blades is to assume that only a single mode is required in 
each direction of bending and further that the mode shapes employed are 
those of a non rotating uniform cantilever beam. 

1 

Another approximate treatment has been suggested in the literature 
which assumes that the blade and hub structural properties can be 
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represented by torsional springs located at the root of the blade. 

The blade is assumed to be rigid and the hub stiffness is represented 
by a set of springs which remain at a fixed orientation to the rotor 
shaft. The blade stiffness is represented by a set of springs which 
rotate as the blade pitch angle is changed. This lumped model is par- 
ticularly convenient for obtaining insight into the manner in which 
various structural parameters enter into the stability analysis of 
flap-lag motion. 

A basic parameter in the flap-lag stability problem is the blade 
pitch angle (or thrust level) of the rotor. Typically, instabilities 
occur as the blade pitch angle is increased^. The importance of the 
blade pitch angle to the instability, as well as the sensitivity of 
the analysis noted earlier, raises the question of whether the two < 
approximate models described above result in a similar functional 
dependence of the structural properties (the stiffness matrix) of 
the blade/hub system on blade pitch angle. 

Physically it would be expected that there would be a difference 
in these two models for the following reason. If it is assumed that 
the mode shapes of the hub/blade system are independent of the distribu- 
tion of stiffness then as the stiffness is varied, the relative deflections 
between segments remain the same since the mode shape does not change. As 
a result, local equilibrium along the blade is not satisfied as stiffness 
is varied and the assumption of a fixed mode shape is equivalent to assuming 
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that various segments deflect proportional to each other independent of 
stiffness. This is in contrast to the model in which the hub stiffness 
is modelled as a pair of torsion springs and the blade stiffness is 
modelled as a pair of torsion springs. This model, taken with the 
assumption that the hub has no mass, determines a relative deflection 
between the hub and the blade which is a function of the relative stiff- 
ness of the hub spring and the blade spring. This spring model is 
referred to as a series spring model implying, that relative deflection 
between the inner spring (hub) and the outer -spring (blade) is determined 
by the relative spring stiffness . • 

It would be expected that the approach in which it is assumed that 

the mode shapes do not depend on stiffness would be analogous to 

the assumption that the relative deflections of the hub and blade are 

\ 

proportional and that the constant of proportionality does not depend 
on the relative stiffness of the hub and blade. This assumption corre- 
sponds to what is referred to as a parallel spring configuration in the 
following. 

The similarities and differences between these two modelling approaches, 
the spring model of Keference 1 and the fixed mode shape model are developed 
in detail in this report. A spring model is developed which is equivalent 
to the fixed mode analysis in the sense that it produces an identical de- 
pendence of the stiffness matrix on blade pitch angle as the fixed mode 
analysis. The pitch dependence of this model, referred to as a parallel 
spring model, is compared to the pitch dependence obtained 'from the series 
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spring model and the differences discussed. 

Also root locus techniques are developed to study the influence 
of various parameters on the stability of the flap -lag motion with 
particular reference to the importance of the difference in pitch 
angle dependence obtained from the two model assumptions. 



EOOT LOCUS STUDY 


Owing to the sensitivity of the problem of predicting flap-lag 
stability of a hingeless rotor blade, an examination of the problem 
employing root locus techniques is highly desirable. This technique 
can be conveniently used if it is assumed that the blade pitch angle 
is small. The following analysis considers, in particular, the sensi- 
tivity of the movement of the flap-lag roots with increasing blade 
pitch to various assumptions which may be empolyed in studying the 
problem. The difference between the results obtained if a flexible 
blade model is used in contrast to a rigid blade with springs at the 
root is considered. The influence of a- simple approximation to the in- 
flow is also examined. Further, the results obtained using a series 
spring model are contrasted with those obtained with a parallel spring 
model, it is shown elsewhere in this report that the parallel spring 
model is .analogous to the results obtained from a flexible blade model 
in which the mode shapes are taken to be independent of pitch angle . 

First the' characteristic equation describing the flap-lag motion is 
developed in a form suitable for root locus studies. 

The perturbation equations for the flap-lag motion given by Ormiston/ 
X 2 

Hodge s/Bousman’ are 


where 


s® + F*s + F- -sF* + F^ 

p P t t 


AP I 

-sC» + C- s® + Cl s + 

P P t t 


Ac| 

- 


/ 


Fp = 1 + i [u)| + H sin® 0] 


[u)J - N sin® e] 


CD 


5 





sin 2 9 


C 






A = 1 + NL sin® 9 


where N ^ R(od^ ~ co®) 


( 2 ) 


L = 



W 





R - as defined in Equation (9) of Reference 2j see Equation (23) 
et.seq. herein. 

The notation R and L are introduced so that the matched stiffness case 

i 

(cOp® = U3g® ) can be conveniently discussed. In addition these definitions 

serve to distinguish between the series spring model as developed by 
1 2 

Ormiston ’ ^and the parallel spring model as developed in a later section 
of this report. It is shown that the parallel spring model is equivalent 
to the formulation with A = 1; i.e,, setting L = 0 in the resulting 
characteristic equation gives the characteristic equation for the parallel 
spring model. 

The remaining terms in the equations of motion are 


F* = T1 (29 - 0^) -2P^ (3) 

= -T| (e - 20^) +2P^ 

C 

= T1 (2 + 00. ) 

Y 

where 'll - g and structural damping has been neglected. 

The coning angle is given by the steady state equations as 
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( 4 ) 





1 

(e - 0i) 



j 

II 

c 


1 

lOJ 

o 



- — - 00. + 1 0.3 
a ^1 8 1 

, 


Therefore the coning angle is equal to 




(9 -0^) -HF^ 

''c 


C5) 


It is assumed that the pitch angle is small such that sin 0=9, cos 0=1. 
Further, it is assumed that the inflow angle at three quarters radius, 0^, 
can be approximated by the following expression 

- 6 I (6) 

noted in an earlier paper by Ormiston^ with 8=1. The constant 6 is re- 
tained to permit examination of this approximation. It can be seen in the 
paper by Qrmiston^ that for a solidity in the neighborhood of .15, 6=1 
is quite a reasonable approximation; however, as the solidity is reduced, 
comparison of this result with momentum theory indicates that a smaller 
value of 6 should be chosen to approximate inflow dependence on pitch angle. 
A 6 of 0.5 will be examined to approximate the low solidity case. 

Further, consistent with the approximation that the pitch angle is 
small, the following approximation is made 


1 

A 


1 

1 + NL sin^ 0 


= 1 


- WL0® 


(7) 


This is a consistent approximation when the pitch angle is assumed small, 
except possibly in the ease of a high chordwise stiffness. 
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The form of the result for the characteristic equation is such that 
it is necessary only to retain the linear term for the dependence of coning 
angle on blade pitch to obtain the dependence of the flap-lag roots on 
pitch angle. Therefore, introducing the notation p^ = 1 + 


^o - ^ - 2^ 


( 8 ) 


The small term depending on blade profile drag has been neglected. 

The various terms in the equations of motion are therefore , approximately, 
using (8) in (2). and (3) 


Fp = P® + W (1 - rxop® ) 0® 

- IT (1 + ) .0® 

Cq = = N 0 

P C 

Fp - n . (9) 

= 2 ne {(1 - (1 - |)] 

Cp = 2Tie + ^3 (1 - |)] 

c; = T1 (2 -^ + I e®) 

q a c: 

The characteristic equation is from (l) 

(s^ + F^s + Fp) (s® + C's + C^) - (sF^ (sC^ “ Cp ) = 0 (lO) 


Using the approximations from (9) and retaining only terms of order 0^ 
for the pitch angle dependence, the following is obtained from (lO) 


8 



(s^ + T) s + p.®) (s® 4- Tl (2 — ) s + w/ ) 

a b 


+ 0' 


jf- s j s= + U (l - i (1 - |)= i (1 - « ) (1 - 6 )) 


+ N 


(n I s . ^ - l) 


-NLtDj.= ((1 + S= + HS + (1 + aOg=))j = 0 


))’ 


(11) 


where the N, L notation has been retained at this point so that the parallel 

and series spring models may be clearly distinguished. 

For the parallel spring model 

L = 0 (12) 

1 

For the series spring model of Ormiston 

N 


L = 


U) .® 


[1 — R ] 

R 


(13) 


The parameter W = R (u)^® - u)p® ) and was introduced to make it convenient 
to examine the matched stiffness case. For the flexible blade case, the 
characteristic equation will of course be of similar form with different 
coefficients as will be noted later. 

The lag damping term ensuing from the profile drag of the blade 
2C^o 

('ll ) has been neglected when it appears added to the flap damping 'll . 

cl 

Now, consider various eases of interest in flap-lag stability. First 
it may be noted that in either of two particular values of the elastic 
coupling paramter R; R = 0 (the principal axes of the blade/hub system 
remain parallel and perpendicular to the rotor shaft and do not rotate 
with pitch) and R = 1 (the principal axes of the blade/hub system rotate 
an equal amount as the blade pitch angle ) , the parallel and series spring 
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models yield the same characteristic equation since L = 0 when R = 1 
for the series spring model. These two limiting cases are considered 
first. 

a.) Ro Elastic Coupling (R = O) 

In this' case R = 0 and the characteristic equation reduces to 

C, 

(s^ + 71 s + p®) (s^ + 71 ( ) s + ) 

+ a|e^ rs{s=-+Tl (l (1 - |f( ^ (1 . |) (1 . 6)) s . 

+ =.o (R = 0) (14) 

The poles of this locus are the uncoupled flap and lag motion. There are 
three zeros, one at the origin, and a complex pair with imaginary parts 
approximately located at the uncoupled rotating flap frequency, . and real parts 
dependent upon the uncoupled rotating flap frequency (p^) and the assumption 
regarding inflow (the parameter 6 ) and the lock number. 

Select the following physical parameters 
P^ = 1.33 

71 = .625 (y = 5) 

’a = 5.73 

= .01 

40 , - . ^ 

such that the results may be directly compared with Ormiston/Sousman. The 

influence of various values of the lag freqqency are of interest so it is 

not selected at this point. With these values (l4) becomes 

(s^ + .625 + 1 . 333 ) (s^ + .002 s + 0)^^ ) 
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+ .3125 6 02 + -^25 (l - ^ (1 - |f + ^ (1 - |)(1 - 8)) s 

+ 1.333|J = 0 (15) 

If 6 = 1 , the quadratic factor for the zeros is equal to 

ts® - .0785 + 1.333] 


If 5 = -. 5 j this quadratic factor is equal to 

[s^ - .352 s® + 1.333] 

Note the strong influence of the inflow model on the location of the real 
part of the complex zeros. Note also that the parameter 6 appears in the 
root locus gain. Denote 
D = T 1 !* (1 - 6)2 (P° - 1 ) -g) 

6 ^ p* 

such that the characteristic equation is ~~ 

• (s® + .6255 + 1.333) (s^ + .002 s + ) 

+ .3125 6 0 ^ [s (s^ + Ds + 1.333)] = 0 (16) 

As an aside, before discussing the root locus>a simple result can be obtained 
for the coupled roots in the special case where " 1 * 333 . Neglecting 

the blade profile drag term, this case results in the characteristic equation 
of the form 

(s® + D„s + p®) (s® + D,.s + p®) = 0 (17) 

X ij 

where 

D (1 + ^) = Dj + Dj, 

Tl I 0 ^ D = D^D^ 
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and since the lag damping is small 

Dj. = T1 ( 1 + ^) 


Dj^=|d9 = 


Showing that a negative value of D leads to an instability in the lag motion. 
For the particular values above the lag damping as a function of pitch angle 


is given by 


. 039 

= - . 088 e® 


(6 = 1 ) 
(6 = .5) 


in the special case where txi^ 


Note the sensitivity of the result to 


the selection of the inflow parameter 6. 

Now to continue with the root locus considerations . Figure 1 shows the 

location of the zeros and poles for = 1.1 and the locus of lag roots. The 

lower figure with 6 = 0.5, corresponding to a low solidity agrees well with 

1 

Figure 6 of Ormiston, et. al. , as would be expected. For simplicity the flap 
mode trend is not shown. The form of this locus will be quite similar for other 
cases since the complex zeros are located at the flap frequency. It can be clearly 
seen that the most critical case as far as destabilizing the flap mode is 
concerned occurs when is in the vicinity of p®. Figure 2 shows the 
influence of lag frequency on the locus. 6 is taken as 0.5 to agree with 
Ormiston 's results presented in Figure 6 of his paper?" 

Note particularly the importance of the inflow approximation (or 
essentially solidity); a value of 6 =1.0 indicates that small pitch angles 
will stabilize the lag mode for = 1.4 where a 6 = 0.5 indicates that 
increasing pitcn angle destabilises tbe lag mode. At = .7 tne trend 
is less sensitive to 6. The matched stiffness case presents a special 


situation discussed later. 
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Recall that there is no influence of the parallel or series spring 
approximation . 

b.) Full Elastic Coupling (R = l) 

Now consider the other limit in which R = 1, The characteristic 

equation in this case is from (11 ), recalling the definition of D, 

2C, 

(s® + T1 s + p®) (s® + T| s + u) ® ) 

+ e® s (s® + Ds+ p®) + ((o| - w®) (T1 I s - l)j =0 (19) 

The zero Ideation for the root locus now depends upon the difference 

between the nonrotating flap and lag frequencies. A locus of zeros with 

variation in o)^ can be conveniently sketched as shown in Figure 3* The 

matched stiffness case is the same zero location as R = 0 for the special 

case when R = R-, and this locus is the same for either the series or 
P t 

parallel spring model. The matched stiffness case is treated in detail later. 
There is a difference between the parallel spring model and the series spring 
model when R* 7 ^ R-. 

P C 

With the zero locations given, the root locus for increasing pitch angle 
for various values of lag frequency can be conveniently sketched as shown 
in Figure 4. Generally in this case where R = 1, the lag mode is stabilized 
by increasing pitch angle. At higher lag frequencies there appears to be a 
tendency towards instability but care must be taken with the approximations 
as noted. 

These two limiting cases of R = 0 and R = 1 serve as a valuable view 
of the limits encountered in the more general case of values of R between 
zero and one. The zeros of the loci- for increasing pitch angle at various 
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intermediate values of R will move from the location shown in Figure 1 
to the location shown in Figure 3* However, it should also he recalled 
that in the cases where a value of R between zero and one is considered 
there will be a difference between the results depending upon whether the 
parallel spring model or the series spring model is employed. Recall also 
that the parallel spring model is analogous to a flexible blade model with 
an assumption of mode shapes independent of pitch. 

c.) Intermediate Elastic Coupling (0<R<1)-. 

In this case the characteristic equation will differ in the 
parallel spring case (L = O) from the series spring case 

L = L (1^) 

The characteristic equation for the parallel spring model is from (11 ) 

2C 

(s® + T1 s + p®) (s® + T) ( — s + U) ® ) 

a Q 

+ 0® |s (s®+Ds + p®) 

i R ^ ^ _ ^^,^8 ) (1 _ R) _ ijjj = 0 (20) 

Note'lhat the zero location bears a certain similarity to the R = 1 
case with the exception of the term 

i (“ c " - V ) (^-^'1 

Without this term, the locus of zeros with R varying would be identical to 
the locus shown in Figure 3 for varying with R = 1. 
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For the series spring model, the characteristic equation is 

C, 

(s® + 11 s + p®) (s® + 11(2 — ) s + U) ® V 

» t / 

+ e® [ ~ s (s® + D s + p®) 


+ R (u)^® - U)p® ) s + (O)^® - U)g® ) (1 - R) - 

( 21 ) 



where there is an additional quadratic factor in the expression for the 

zeros, the underlined terms, which arise in the series spring model, but 

are not. present in the parallel spring model. 

The series spring model will yield results which are the same as those 

1 

obtained by Ormiston and shown in Figure 6 of his paper. The parallel spring 
model will yield quite different results. 


First consider the case in which the lag frequency u)^ = 1.1. Substituting 


this and the other values 


P® = 1.333 a = 5.73 D = -.078 6=1 

11 = .625 = 0.01 D = -.352 6 = .5 


the characteristic equation becomes 


(s® + .6255 + 1.333) (s® + .002 s + 1.21) 

+ .3125 6 6® j^s (s® + D s + 1.333) 

+ .88r (s + ^ (.88(1 - R) - 1)) 


-(1 - R) (8.4l6) (1.275 s® + .625 s + 1.666) j = 0 

6 J 
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where the underlined terms are dropped for the parallel spring model. In 
this example we select 6=1, giving the following characteristic equation 


(s^ + .625 s + 1.333) (s® + .002 s + 1.21) 

+ .3125 0^ |^(s) (s^ - .078 s + 1.333) 

+ .88 r |s + 3.2 (.88(1 - R) -1)| 

- 10.73 (1 - R) (s® + .490s + 1.306)j 


= 0 


Now then to compare the series and parallel models the zero locations for 
the series and parallel cases are shown in Figure 5. For reference, the 
locus of zeros given in Figure 3 = is included as this represents the 

limiting points. 

Note the sensitivity of the zero locations to the use of the series or 
parallel model . 

Figure 6 shows the equivalent diagram for = 1.4 where the character- 
istic equation is for 6=1, 

(s^ + .625 s + 1.333) (s^ + .002 s+ 1.21) 

+ .3125 0® (s) (s® - .078s + 1.333) 

+ 1.63 R |s + 3.2 (1.63 (1 - R) - 1)| 

- 18.43 (i - R) (s® H- .535 s + 1.426)J = 0 

Figure 6 shows the locus of zeros for = 1.4. It can again be seen 
that the zero location is quite sensitive to the model employed. Figure 7 
then shows the difference in the departure angles for these two models for 
0)^ = 1.1 and u) = 1.4. The limiting eases (E =0, R = l) are the same. 
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however, the departure angles differ markedly again indicating the 
sensitivity of the dependence of the lag modes to modelling assumptions 


Recall also that particularly at = 1 . 4 , there is a sensitivity to the 
inflow parameter 6 . 

The parallel spring model also indicates one interesting special case 
corresponding to a particular value of R, in which the locus configuration 
is similar to the R = 0 case. Note that this result will not agree with 
the series case. The particular value of R is 


(u)^® - ) (1 - R) = 1 

or 



For this specific value of R, the parallel spring model gives zeros at 
s (s® + D' s + ) =0 

Thus this special case will always give a locus which is similar 
in shape to that obtained for R = 0. The complex zeros are located at 
the lag frequency and in the right half plane if D is negative, 
d. ) Matched Stiffness = ixi^ ) 

In the matched stiffness case the parameter N, introduced 
earlier is employed, rather than R which becomes infinite except in the 
case R^ = Rq . Therefore 

C P 

N = R(u)^® - 

If Rg = Rq, N = 0 and the matched stiffness case will be the same as R = 0 
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If Rg ^ Rg as would be more typical of the rotor blade, then there will 

be a difference in the locus from the R = 0 case. The characteristic 

equation in this case is to ® = (Oq® = u;P 

£ P 

C, 

(s® + 71 s + p®) (s® + T| (2 ~) s + to®) 



where again the underlined terms appear in the series spring model and 
not in the parallel model. 

Since in this case 




AR = Rp - R^ 


and a typical rotor blade would have a considerably higher chordwise 
frequency than flapwise, positive values of N are of interest. 

The departure angles are shown in Figure 8 for the matched stiffness 
case again using the previous parameters 

Y =5 .01 

0)^ = .333 a = 5.73 

for W = 0, where the blade stiffness is matched and the results are the 
same as for R = 0, and for AR = .75 which would be more typical of a 
rotor blade. 

The coupling in this case has a strong stabilizing tendency on the 
flap-lag motion. Again the two models exhibit a difference. AR could 
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be somewhat larger than shown. It is eq^ual to (in terns of frequency) 



and so is essentially limited to a maximum value of about one. Negative 
values of AR would be highly unlikely physcially as this would imply 

that CO- > (0^ . 

Pb ■ 

This completes the discussion of the rigid blade models and the 
comparison of the series and parallel spring models . Now the rigid and 
the flexible blade models will be compared . The flexible blade model 
considered here (See Appendix I for a discussion of the flexible blade 
equations of motion. ) assumes that the mode shapes are independent of the 
blade/hub mass and stiffness distribution and as shall be seen is structurally 
equivalent to the parallel spring model discussed above. 

If the flexible blade model is used to determine the stiffness terms 
in the equations of motion and consequently the equivalent terms in the 
rigid blade model, then these terms for the flexible blade model will be 
identical to those of the parallel spring model. 

The differences in the equations of motion in the flexible blade 
case as ccanpared to the rigid blade case, will arise therefore in the 
aerodynamic terns and the coning angle . The coefficients of these terms 
will be changed from the rigid blade values depending upon the mode shape 
employed. 

In terms of the previous discussion the following terms in the 
characteristic equation are unchanged therefore. 
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Fp = p® + N sin^ 6 
- W sin^ 0 
= Cp = I sin 2 6 

where A is jtaken as 1 since mode shapes independent of pitch angle are 
employed* p®and q® are the rotating flap and lag frequencies and include 
the Southwell coefficient which no longer equals 1 as it does in the rigid 
blade ease. 

The remaining terms are as follows. A uniform mass blade is considered. 
RIGID BLADE FLEXIBLE BLADE 

= 11 T| a^ 

F • = 11 (20 - 0^) - 2 ^ ^ 

Cr = - Tl (0 - 20. ) + 2 P = - 11 (a^^ + 21 a^) + C w° 

p 1 o p o 

c; = T| (2 — ) + 00 . c; = 11 (2 — a^^ - X 

B. X Q Qj 

Table I gives the values of the various mode shape parameters given above 
when a Duncan polynomial is used for the first mode shape. 

•Substituting numerical values for an untwisted rotor blade the following 
results are obtained. The terms for the rigid blade can be obtained directly 
from the mode shape integrals with ijt ~ x. See Appendix I. 


RIGID BLADE 

FLEXIBLE BLADE 

F^ 

1.070 11 

F^ = H (20 - 0^) - 2 p^ 

11 [2.14O0 - 1.335 0t] - 2.110 w° 

= -H (0 - 20^) + 2 p^ 

- T| [1.070 0^ - 2.6.70 0t] + 2.110 w° 

C 

= T1 [2 -|^ + 00^] 

G 

n [(1.070) (2 -p) + 1.335 90t] 

3. 


20 



The coning angle is given by 


\ (e - 

p 


-o ^ U (.Xrtw ^ ^ 
° p= - 


Pc = ^3 


^ [1.1688 e - i.50i 4 0f] 


if 

where 0 t = -^ 0^^^ in the flexible blade equations, giving for the flexible 
blade the following coefficients 


Fp = 1.070 T1 

F* = T1 [2.i 40 0 - 1.001 0.] - 2.110 w° 

Cp = - T] [1.070 0 - 2.003 0j^l + 2.110 w° 

C 

C; = 11 [1,070 (2 —) + 1.001 0 0.] 
t a 1 

w° = 3. [1.1688 0 - 1.126 0.] 

p2 

6 

Introducing the approximation that 0^ = where here for simplicity the 
6=1 only case is examined. 

f; = (i. 64 o - -•— ) lie 

^ p" 

= (-.068 + lie 

P P^ 

F* = 1.070 11 

C 

c: = 11 [1,070 (2 -^) + .500 0®] 

-O ^ .6058 119 
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For comparison purposes, ’the rigid blade terms are 


C, 

G' [2-~ + .5 e^] 

Fi = 11 9 (1.5 - — ) 

= H 0 (i) 

.5H e 

3 = 

0.2 


The significant change appears through the coning angle affecting 
primarily the coupling terms F* and . 


(s^ + F^s + Fp) (s® + C^s + C^) - (Cp -Cp) (F^s - F^) = 0 . (lO) 


The R = 0 case will be examined for these two models when R = 0, F- = Cq = 0 

’ C 3 

and the pitch dependence appears in the term, the lag damping, and the 
coupling term -C^ F^ s^. 

The characteristic equation is therefore, with the 9 dependent terms 
separated out equal to 

C, 

(s® + F^s + Fp) (s® + 11 (2 ~~) s + C^) 


+ 0® I s (s^ + [F^ - ^ ® ° 

® 2 

The term which determines the real part of the complex zeros is essentially 
the only place in which the difference in the two models appears. This is 
the term which was denoted as D earlier. 
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D = F^ - 


!il 

9=1 


For the rigid blade model 


D = T1 - ^ (1.5 - — ) = T] ~ ~ 

This shows clearly the limits in flap frequency found by Ormiston^ If 
1 < p® < 2 then the complex zeros for the pitch angle locus are located 
in the right half plane and for rotating lag frequencies in the vicinity 
of the flap frequency, the lag poles will be attracted to the right half 
plane by these complex zeros leading to instability with increasing pitch 
angle . 

For the flexible blade (with Duncan polynomials for mode shapes) 


D = 1.070 7] - 2T1 (1.640 - _ .068) 

p® p® 

D = 1.293 71 - 1>119) 

Thus the flexible blade model with Duncan polynomials as mode shapes will 
give a somewhat different frequency range for instability. Instability 
would be expected for 1.12 < p® 2.26, 

The location of the complex zeros as a function of p for the rigid 
and flexible models is shown in Figure 9* 

The inflow model will affect these results. If 6 =0.5 then the 
rigid blade model gives for D 

p . 4.5 H -2) 

,,4 
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It is interesting to note that selecting 6 = 0.5 does not change 
the range of flap frequency over which an instability would be expected 
but it does increase the movement of the complex zeros into the right 
half plane. 

In the flexible blade case with 5 = 0.5> 

D = 5 269 T1 — 7. .^*35) (p — 

There is again an increase in the variation of the real part of the com- 
plex zeros with changes in flap frequency and the frequency range over 
which the complex zeros are located in the right half plane is about the 
same as in the case where 6 =1. Recall that 6 also appears in the root 
locus gain so that these changes in the zero locations do not directly 
imply a proportional decrease in the stability of the lag mode. 

In the particular case here note that for flap frequencies less 
than\/T^12, in the case where 6=1, the flexible blade model will not 
predict an instability with increasing pitch angle while the rigid blade 
model will owing to the fact that the zeros for the flexible blade will 
be located in the left half plane. 

Thus the nature of the stability of this system appears particularly 
sensitive to modelling assumptions for the elastically uncoupled case, 
and the results will depend upon the assimiptions employed. Physically 
this can be seen to arise from the fact that the real part of the complex 
zeros is located by a term which is the difference of the product of the 
flap damping and the part of the lag damping dependent upon pitch and 
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the rate coupling terms arising from aerodynamics and coriolus terms, i.e., 


The second term subtracts from the first and so this term may- 

be of either sign and is quite sensitive to the assumptions employed in 
studying flap-lag stability. This is particularly the case when the un- 
coupled lag frequency is in the vicinity of the flap frequency such that 
the lag roots are directly attracted to these zeros. The Important role 
played by the dependence of the lag damping on blade pitch can also be 
noted from this expression. If 6 is set equal to zero, the lag damping 
becomes independent of pitch angle and there are two zeros at the origin. 

The root locus angle condition will be 0° for 1 < p^ < 2 and l80° for p^ 
outside this range, indicating a more severe instability with increasing 
pitch in the range of 1 < p^ < 2. The .presence of the lag damping" depen- 
dence on pitch acts to limit the instability with increasing pitch to some 
maximum unstable damping as shown by the zero locations in Figure 9. 

It is ^pected that the R = 1 case is not so sensitive as the zeros 
in this case are moved into the left hand plane by the additional structural 
terms introduced by elastic coupling as can be seen from Figure 3. 

For small values of R it would be expected that the sensitivity to 
modelling assumptions will exist. It has already been shown that for small 
values of R, considerable differences exist between the use of a parallel 
spring model, corresponding to a fixed mode shape assumption and the series 
spring model proposed by Ormiston which physically appears to correspond to 
a more refined assumption regarding the variation of the mode shanes with 
blade pitch. 
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To summarize this discussion it can be seen that for a hingeless 
rotor blade with no elastic coupling the results of a flap-lag stability- 
investigation are sensitive to the details of the assumptions made in the 
physical model. 

For small values of the elastic coupling parameter a similar sensi- 
tivity exists which depends upon the nature of the structural model as 
well as the assumption regarding the inflow. 

The fully elastically coupled blade (R = l) does not appear to be 
particularly sensitive to modelling assianptions on the basis of these 
investigations and the same results are obtained from either the series 
spring model or the parallel spring model. 

Intermediate values of R appear to give rise to a situation in which 
the results of a stability investigation are very sensitive to the modelling 
assumptions. There appears to be less sensitivity for blades which are soft 
in the chordwise direction than for blades which are stiff in the chordwise 
direction judging by comparison of the trends shown for a matched stiffness 
blade compared to the stiff inplane case. As the chordwise frequency is 
increased there appears to be an increasing sensitivity to the use of a 
series spring model as compared to the parallel spring model as can be seen 
by comparing Figures 5 and 6, implying that for intermediate values of R, 
the details of the modal analysis for a flexible blade are significant. 

The root locus technique appears quite valuable in pointing out some 
of the sensitivities involved in this problem. It would be expected from- 
these results that while the general shape of a stability diagram such as 
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shown in -Figure ^ of Ormiston's paper^ illustrates the nature of the 
stability boundaries for flap-lag stability, the precise values will 
be sensitive to the details of the model. Further it would appear that 
for the results shown in Figure 6 of Ormiston's paper, the trends shown 
for 0 < E < 0.4 are also very sensitive to the particular analsrbical 
model used. 

To make this point more concretely, Figures 10 - 12 have been 
prepared to show the locus of roots of the lag mode as a function of 
blade pitch for three values of the uncoupled lag frequency, based on 
the parallel spring model. These figures correspond to the root locus 
sketches in Fig^^res 7 and 8. The figures may be directly compared with 
Figure- 6 of the paper by Ormiston and Hodges and illustrate the 
differences between the parallel spring model and the series spring model. 

The two models give quite similar results for = 0.7 and considerably 
different results for = 1.1 and l.U, Note the decreased sensitivity 
to pitch angle predicted by the parallel spring model as compared to the 
series spring model by Ormiston and Hodges. 
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THE PARAMETER R 


The parameter R was introduced by Ormiston to represent the elastic 
coupling of the blade/hub system. It can be conveniently interpreted as 
being directly related to the rotation of the principal axes of the blade/ 
hub system. It can be shown that the rotation of the principal axes of 
the blade/hub system Yp is related to the blade pitch angle and the 
parameter R by the following equation (see pp. 30 and 4l) 


tan 2y 

P 


R sin 2 9 

r""cos" 2 0 + (i '-r) 


( 23 ) 


For small pitch anges this reduces approximately to 



Figure 13 shows a graph of -the principal axis inclination as a function of 

blade pitch angle for various values of R. Note that when R = 1, Y = ® 

6 

for all values of 0, and when R = .5? Yp= 2’ 

'In this investigation the range of R examined was between 0 and 1. 

By placing R in a somewhat different form, an estimate of the possible range 
of E can be seen. As defined by Ormiston 


where 


R = 





"p 




K 


C 




s 
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where cOp ® and U)^ ® are the natural frequencies of the blade rigidly 
mounted. R can be expressed as 



It can be seen that R = 1 if the hub mounting has no influence on the natural 
frequencies such that 

-"c 

R = 0 if the blade itself is of matched stiffness, which is a physically un- 
likely situation unless the blade is softened at the root in the chordwise 
direction outboard of the pitch change bearing. R is also equal to zero 

if the blade is articulated in the flapwise direction such that uJo = 0. Thus 

P 

in the root locus studies described in another section, the case R = 0 implies 
that outboard of the blade pitch bearing, the nonrotating frequencies in flap 
and lag are equal. 

The case R = 1 can be achieved by a number of combinations including 
the situation in which the hub is very stiff such that the blade/hub frequencies 
are the same as the blade rigidly mounted. If the blade is not softened 
outboard of the pitch bearing then it would be expected that (u„ /m- would 
be small leading to a value of R which depends upon the blade/hub system 
frequency ratios and the reduction in the flapwise frequency ratio produced 
by the hub. This gives the possibility of R values larger than 1 particularly - 
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for a blade/hub system which is soft in the chordwise direction. 

Negative values of R, although theoretically possible from the'fom 
of the expression, appear to- be highly unlikely in practice. 
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COMPARISON OF VARIOUS STIFFIjESS MODELS 


In this section, the form of the stiffness terms in the elastic 
blade equations are examined for the case in which it is assumed that 
the bending of the blade in the flap and lag directions is represented 
by a single mode and also that the mode shape is independent of the 
stiffness distribution between the blade and the hub. Of particular 
interest is functional dependence of the terms in the stiffness matrix 
on blade pitch angle. In addition two other simple models of the blade/ 
hub system are examined. The hub and blade are modelled as rigid members 
with the hub stiffness represented by a spring at the root and the blade 
stiffness by a spring between the hub and the blade. Two alternate 
assumptions examined; one in which the blade deflection is assumed to 
be proportional to the hub deflection and independent of stiffness of the 
spring, and the second in which the relative deflection between the blade 
and the hub is determined from the condition of equilibrium of the hub 
segment^. These two models are referred to as the series spring model 
and the parallel spring model as well will be explained. The blade 
pitch dependence of the stiffness matrix for these two models is examined. 
It is shown that the parallel spring model yields a pitch dependence of 
the stiffness matrix identical to the modal analysis with fixed mode 
shapes and the series spring model yields different results, in general, 
as to the manner in which the stiffness matrix varies with blade pitch. 

First consider the elastic stiffness terms as given by the flexible 
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blade equations of motion presented in Appendix I. The terms of Interest 
here are: 


r.^w0 

„VW0~ 



K 

K 


V 

^■vwQ 

,,WW0 



K 

K 


w 

w 




= J <3^+1 (®V ■ ® (’■'v ^ 

o o ^ 

K™® = r El /(ilr//)® dr + (El , - El ,) sin® 0 dr 

o y ^ *^o ^ ^ 

R 

^vwe = [EI^, - Ely,] SiB^ (,,» ar (£4) 

Fow , if it is assumed that the mode shapes are independent of blade 
pitch, that is, they are the mode shapes calculated at zero pitch then 
the functional form of these stiffness coefficients is as follows 


f El , (t;')® dr - sin® 0 f (El , - El / ) (^'')® dr 
o o 

^ww0 ^ j El /(t];'^)® dr + sin® '0 f (EI / - El / ) (ilf ^ ' )® dr • 
fs y ^ r, ^ y ^ 


^ R 

= sin 0 cos 0 r (El / - El , ) (tfr ' ' ijf ' ' ) dr 
J ^ ^ z y ' ' 

If the rotor blade of interest is considered to consist of two 
portions, an inner segment denoted by the subscript h which does not 
rotate, extending to a radius r^^ and an outer portion denoted by the 
subscript b which rotates with the pitch angle 0 , the terms are 
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= ; El , (t'')= dr + f El dr - 


‘h 


R 


z'h 


2 b 


'll 


R 

- sin® 0 r (El / - El /, ) (t'O^ ^ 


K™®. J El ,,(♦;')= dr. J El. (♦;') 


h 


R 


y 


dr 


'h 


R 


+ sin® 0 r (El /, - El /, ) dr 

Jrj^ ^ z b y 


K 


•VW0 


R 


= sin e cos 0 r (El /, - El /^) (ilf'' A^') d.r 

z b y b^ ' 


•h 


(25) 


If it is further assumed that the mode shapes are the same for both 
directions of bending such that 


\jf = ijr =1], 


and using the following notation 




These coefficient beccane 


= k + k , - sin® 0 (k , - k . ) 
zh zb ''zb yb^ 


,ww0 




zb yb ‘ 


(26) 


j^w0 _ sin 0 cos 0 
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■where the stiffness coefficients, the k’s are independent of pitch angle. 
Note that the effective spring constants of the hub and k^^ appear 
as a sum. This result is of a form that would be expected for two 
springs with spring constants k^j^ and k^^ connected in parallel. 

This result can be -written in matrix notation by defining three 
matrices 


k„ = 



0 


k 


■yh 



0 


cos 9 sin 0 


|-sin 9 cos 0 J 

The stiffness matrix may be written as 

[k] = [k„] + [0]-^ [k^l [0] 


(27) 


where all of the effects of blade pitch are contained in the rotation matrix 
[0]. Thus this displays in convenient form the stiffness matrix for the 
flexible blade model in the case in which it is assumed that the mode shapes 
are independent of pitch and are the same for both deflections. There are 
five parameters involved in this form to define the coefficients in the 
stiffness matrix. The four stiffness coefficients and the pitch angle 0. 

The number of independent parameters in the stiffness- matrix can be 
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reduced by one by introducing the concept of a principal axis, that 
is by finding the transformation which diagonalizes the above matrix 
and then expressing the rotation of the principal axis in terms of the 
pitch angle. Consider therefore, a rotation matrix [y] which diagonalizes 
k by a similarity transformation ■ 


[K] ^ [y]-" Ck] [y] 


If the matrix [k] is denoted as 


k H 


•11 


■13 


k 


13 


■33 


the rotation which diagonalizes the matrix k is y where (cf. ( 23 )) 


2 k 


tan 2y = 


13 


P kii - k 


33 


Ak, sin 2 0 
b 

Ak, + Ak, cos 20 
h b 


(28) 


where 




it E (k ^ - k 
h zh yh^ 

y is the rotation of the principal axis which is given above as a 

XT 

function of the blade pitch angle and the various stiffnesses. It can 
be noted that if the blade stiffness is equal in both directions then 
the principal axis does not rotate with blade pitch. If the hub stiff- 
ness is equal in both directions then the rotation of the principal 
axis is equal to change in pitch angle. 
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The stiffness coefficients in the diagonalized matrix are: 




K 


= I (s + E V + 1 


■) 


( 29 ) 


K 


- = 2 (S + S is) - I ( 


^ Ak^ + Ak^ cos 20 


cos 2y 


•where 


*lyn) 


^■s “ (’'zb ■'■ Vb> 


It can be seen that if the blade stiffnesses are equal then these 
two stiffness coefficients are independent of pitch angle. Similarly 
if the hub stiffnesses are equal these coefficients are independent of 
pitch. That is, in both these limiting cases, the two nonrotating 
natural frequencies of the blades will be independent of blade pitch. 

Now these results are compared to rigid blade-models with root 
springs which. may be assumed to represent the structural characteristics 
of the blades for use in an analysis of flap-lag stability in a somewhat 
simpler fashion. It will be shown that the spring model prposed in 
Eeference 1 gives rise to a different pitch dependence from that given 
by the stiffness matrix of equation ( 27 ). 

Therefore, two spring models will be considered, the one proposed by 
Ormiston (Reference l) which will be referred to as the series spring model 
and a second one referred to as the parallel spring model. ^ These two models 
will be compared to the results obtained from the flexible blade assuming 
that the mode shapes are independent of pitch angle. Particular interest 
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centers around the manner in which the results depend upon pitch angle. 

Consider the structural properties of the hub to be represented by 
two orthogonal springs which produce moments proportional to deflection 
in an axis system parallel to the perpendicular to the shaft. The struc- 
tural properties of the blade are represented by two orthogonal springs 
outboard of the hub which are aligned parallel and perpendicular to the 
chord line of the blade (precisely speaking the principal axes of the 
blade). In order to make the notation somewhat more compact the sub- 
scripts 1 and 2 refer to the -two directions. A primed quantity is 
referenced to the blade axis system and an unprimed quantity is referenced 
to the hub axis system, Therfore, the stiffness matrices are 



where these stiffness coefficients are independent of pitch.- The de- 
flections about these two axes are denoted as follows 



and of course denotes the hub deflection in the blade axis system and 

{0 b} blade deflection in the hub axis system. These angles are assumed 
small so that it is not necessary to introduceEuler angles. Figure l^J- 
shows the definitions of these quantities. These quantities are related 
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by the rotation matrix 


[0] = 


cos 9 
-sin 0 


so that 


sin 0' 
cos 0 


{ 0 '.] = [e] f 0 s} 

M =t9] {«„] 

The potential energy of this spring system can be written as (Meirovlteh^) 

V = I j'' [K„] + I [Kj] c 0 i - «'„} ( 31 ) 


where the superscript T indicates the transpose of the matrix. 

The coordinates in the second term are expressed in terms of coordinates 
in the hub system using the transfonnation (30) 

V = I {0h}^ [K„] {0h} + I {0b - 0h}’^ K^J [0] {08 - 0h] (32) 

where the notation 

[Kb] = [0]^.[k;] [9] 
is introduced such that 

V = I {0„}^ [K„] {0,} + I {0B - 0h1^ [Kb] {0b - 0h] ■ (33) 

In this expression for potential energy there are four degrees-of-freedcmi 
the blade deflection in two directions {0 b} and the liub deflection in two 
directions {0^}* The two spring models to he discussed involve different 
assumptions as to the relationship between these degrees -of -freedom. Two 
spring models are developed to show clearly the analog between the fixed 
mode analysis and a rigid blade model with root springs. The parallel 
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spring model is developed with the specific objective of obtaining a 
stiffness matrix which varies in the same fashion with pitch angle as 
the fixed mode analysis (Equation (27)). 

The parallel spring model consists of making the assumption that the 
hub deflections are proportional to the blade deflections 

{0„} = [a] {0 b} (34) 

such that 

V = .| {0.]’^ m’’ [k„] [<»] [0b3 + 1 [0.}’^ ([I] - [o])’^ [K,] ([I] - M) {0,} 

(35) 

and the stiffness terms are therefore 

- [c»]® [Kh] to] {0J + ([I] - to])’' [K.] ([I] - w) {0,} (36) 

In the parallel spring model it is further assumed that the a matrix 
involves only one constant such that 

[a] = a [1] . (37) 

and therefore, the stiffness terms for the parallel spring model are 

{||g}p = a® C^h] + (1 - a)® [KbI (38) 

where the bars are introduced to indicate that this is the result for the 
parallel spring model. The pitch angle dependence is contained in [Kg] 
only where 

[Kg] = [0]“^ [K'g] [0] 

Thus the assumption that the hub deflections are proportional to blade 
deflections and are independent of stiffness 3d-elds a stiffness matrix 
that has the identical pitch dependence as the flexible blade model with 


39 



fixed mode shapes (Equation. (27)). 


The series spring model is obtained by allowing {0 g] and {0 to 
be independent and determining the relationship between these two variables 
from the potential energy estpression assuming that the inboard segment, 
i.e,, the hub, has -no mass such that in effect the matrix [a] is found 
- from the following conditions : 


^ 0 „^ 


0 


From (39), using (33) 

- [K„] {0„} 


- [Kb] {0b - 0h} = 0 


(39) 


(40). 


Therefore 

{0h} = [Kh + Kb]“^ [Kb] {0b} (4i) 

that is for the series spring model, cf . (s4) and (4l),, 

[a] = [Kh + Kb]""- [Kg] (42) 

Also 

iaf = [KbI [Kh + Kb]"^ (43) 

since [Kh + Kg] and [Kg] are symmetric. Now from ( 36 ) 

{^^} = {[af [Kh + Kg] [a] + [Kg] - [af [Kg] - [Kg] [a]) {0 b} 
using (42) and (43) 

= [Kg] [Kh + Kg]-^ [Kg] + [Kg] - 2 '[Kg] [Kh + Kg]“^ [Kg] {0g} 

= [Kg] ([I] - [Kh +Kb]-^ [Kg]) {0g} 


4o 
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Hence 


^5(6 ^ ■■ Kg] ^ [Kh] {0b} (^^) 

s 

Thus this is the form of the stiffness terns for the series spring model 
•where again the s-tiffness dependence on pitch angle is given by [Kg] which 
is equal to 

[Kb] = [er^ [k;] [0] 

Thus, to summarize, the stiffness matrix for the parallel spring model is 

[K]^ = a® [K„] + (1 - a)® [Kg] (45) 

and the series spring model 

[K]g = [Kg] [Kh + KgJ-i [K„] (46) 

To show the manner in which the series stiffness matrix varies with pitch 
the series spring model can be ■written as follows, noting that det Kg = 
det Kg since [Kg] and E'Kg] are related by a similarity transformation. 




[K„] + , fKg]] 

^det (K„ + k'b) det (Kh + Kg') ^ 


I ■] , ^Ku AKg .2 a 

II + — ELz_J sin 0 

det (K„ + Kg' ) 

(47) 


It can be seen that the first term in parentheses, (4?), gives the same . 
pitch dependence to the stiffness terms as found in the parallel spring. 


The second term in this e 2 {pression gives rise to an additional pitch de- 
pendence which is a function of the product of the blade and hub stiffness 


differences . 
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New. compare these results to those obtained from the flexible blade 
model with fixed mode shapes which are independent of pitch. This stiff- 
ness matrix-^is of the form 


[k] = [kn] + [0]-" [k;] [6] 


(27) 


The same dependence with pitch angle is found as for the parallel spring 
model as would be expected physically from the assumption of fixed mode 
shapes, see (45). The matrix [Kg] has the form 


K, 


[Kb] = 


11 


- AKg sin® 0 AKb sin 0 cos 0 


AKg sin 6 cos 0 


Kg' + AKg sin 0 

2S 


(48) 


Thus if the blade stiffnesses are equal then' all three of these expressions 
give stiffness matrices which are independent of blade pitch. If the hub 
stiffnesses are equal again all three of these results agree as to the 
pitch dependence of the stiffness terms . In the general case in. which 
both stiffnesses are unequal the series spring model yields a different 
variation of stiffness with pitch from the other two. 

It is further interesting to note that both the series and the 
parallel spring models give the same principal axis dependence upon pitch 
since the difference between the parallel spring model and the series 
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spring model is multiplication of all terms in the series spring stiffness 
matrix by the single factor " 


1 + 


AK„ 

det (K„ + k'b) 


sin^ 0 


The equivalence between the various models is as follows. For the elastic 
model to be equivalent to the parallel spring model, we require 


kn = 0-2 K„ 

11 11 

2S S2 

k'e = (1 - k'b 
11 11 

k's = (1 - cv)2 k'b 

32 22 

where the parameter a is arbitrary.. 

For the elastic model to be equivalent to the series spring model, noting 
that 

det Kb det K„ 

= 6 e and =(l-6) (l-e) 

det (Kh + Kb') det (K„ + Kg) 


where 


k1 


6 = 


11 


Kp' + K, 


11 


11 


We require 


k„ = 6 e Kh 
11 11 


Kp' 


22 


e = 


Kp' + K, 


22 


22 


k|^ — 6 e K^^ 

22 22 


( 50 ) 


k; = (1 - 6 ) (1 - e ) Kb' 


11 


11 


kg' = (1 - 6) (1 - e) Kb' 
22 22 


- h3 



Recall that this equivalence does not account for the multiplying 
factor noted above which is a function of the pitch angle. This 
equivalence will give the same principal axis direction for either 
model. However, the dependence of natural frequency on pitch angle 
will be different. 

These equations give the following relationships between spring 
constants . 


32 

Kh = ) Kh 

32 11 


22 


Kb = (~ ) K„ 


11 


kc 


11 


22 


(51) 


, 22 
Kg = (-— ) K„ 
- 22 k/ 


K, 


11 


( 


Kh +ka +ke % 

11 11 ^ 22 ^•s) 


11 


k ' 


22 


Thus given the results of the flexible blade calculation, equations ('25)- (27), the 

I ! 

spring constants to be used in the) series ispring model are determined by 


K„ = (k„ + ksM 1 + ^ 

11 11 11 \ 


22 


22 


(52) 


kp 


Kh = (k„ + ka 0 1 + 


22 


22 


22 


11 


11 


4U 



ss 


- (^H + kg ^ ) 11 + 

11 11 11 \ 1, ^ 

R 


Kb" = (kH + ks ') I 1 + 


3S 


11 


S2 


82 88 


11 


The correspondence between the notation used here and that of Orraiston/ 


Bousman is 


= (1 - 6 ) 
Ro = (1 - e) 




K's K„ 

11 11 


K/ + K, 


( 53 ) 


11 


11 




K'b. K„ 

88 82 

Kb' + Kh 
88 82 


1 2 

Expressing Ormiston’s parameters ^ in terms of the elastic model with 
fixed mode shapes (equations (25) - (27)) 


^b" 


E. = 


88 


83 28 


, kp 




11 


kg + kn 
11 11 


u)® - (kg + k^ ) 
^ 11 11 

(Up = (kg + k,^ ) 


88 88 
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k/ - ke' 


R = 


11 


32 


(k/ - kg' ) + (kH - k„ ) 


(54) 


11 


33 


11 


23 


The factor missing in the elastic model with fixed mode shapes is 

(kn -kn)(kB -kg) \-l 

11 33 11 22 

1 + sin® 9 

(k^ + kg ) (k^ + kg^ ) 

11 11 23 22 


The factor within the parentheses is A in Ormiston's notation. 

Using the above equivalence^ consider the nonrotating natural 
frequencies and the differences which arise from the two approaches. 
Recall that both models give the same rotation of the principal axes 
with pitch angle. 

The dependence of the nonrotating natural frequencies on pitch 
is given in the series spring case as 


“»1. % <«3= 


-(2(u® 0^®) 


(u)| + ) ± (uJp ^ A-4(r)( 1-R) sin 2 0 

(55) 


In the parallel spring (elastic model with fixed mode shapes-) case? 


cu 


3 




1 - 4 (r) (1-R) sin® e 

(56)- 


These differ by the multiplying factor noted. 


For a helicopter rotor where the blade pitch angle is small, it may be assumed 
that 

4(r) (1-R) sin® 0 « 1 


and the natural frequencies given by the two cases are approximately 
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(1-R) e® 




1 -(^ - 1 ) 
P 


UJi: 


1 “(iU^ - 0 (E) (1-R) e® 

c 


0)^ ®, u) ® = U) 2 

^ P ® P C 


1 +(5%- - i) (R) (1-J^) 6® 

c 


U). 


UJ. 


i + (i - .1) (R) (I-E) 0 


The principal axis rotation given by both approaches is the same using the 
equivalences noted above. In terms of the elastic coupling parameter R, 


tan 2y = 


R sin. 2 e 


p R cos 2 0 + (1-R) 


The frequencies given by each approach in general is different. In the 

/ “’r \ 2 , / 

particular case where the ratio is near to one, that is xf j 


1 + T1 where T| is small compared to one then the same frequency dependence 
with pitch is obtained from both methods. 

Thus using the equivalences given above to determine the corresponding 
parameters each model will give the same principal axis rotation, but the 
parallel spring and elastic blade model with mode shapes independent of pitch 
will give a different frequency variation with pitch from the series model. 

Thus experiments to be described in a later report will not distinguish 
between the two approximations when the principal axis direction is measured. 
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Also if pitch angles typical of a helicopter are examined, only small 
changes in the frequency will be noted and it will be difficult to 
distinguish between these two approximation, series vs. parallel spring 
(or elastic blade model with mode shapes independent of pitch) from 
experimental results. 
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CONCLUSIONS 


As a result of this analytical study of flap-lag stability in 
hovering flight the following conclusions can be drawn: 

1. ) Prediction of flap-lag stability is very sensitive to the 

assumptions used in modelling the aerodynamic and elastic 
characteristics of rotor blades. 

2. ) The equivalent lumped spring/mass model, determined from 

a single mode representation of the elastic deformation 
of a rotor blade/hub system in which the mode shapes are 
assumed to be independent of blade pitch, gives rise to 
a different dependence of the blade/hub stiffness matrix 
on blade pitch angle from the lumped model proposed in 
Eeferences 1 and 2. 

3. ) These two lumped spring/mass models, referred to as the 

parallel spring model (equivalent to the. single mode 
analysis) and the series spring model (References 1 and 
2) give identical results for the variation of the charac- 
teristic roots of flap-lag motion when the elastic coupling 
parameters R = 0 or 1 and different results for other 
values of R. 

4. ) It is possible to determine the equivalent lumped spring 

constants for either the parallel or series spring model 
from the results of an elastic analysis of a rotor blade. 
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5. ) The elastic coupling parameter R can be conveniently 

interpreted as a measure of the rotation of the 
principal axis of the blade/hub system with blade pitch. 

6. ) Experimentally, the elastic coupling parameter R can be 

more precisely determined from principal axis measurements 
than from frequency measurements. 
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FIGURE 1. ROOT LOCUS FOR INCREASING PITCH ANGLE. NO 
ELASTIC COUPLING (R=0), a>f =1.1 . INFLUENCE OF 
INFLOW APPROXIMATION. ^ 
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FIGURE 2 . ROOT LOCUS FOR INCREASING PITCH ANGLE. NO 
ELASTIC. COUPLING (R=0). VARIOUS 
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p2 = 1.33 
/ = 5 


FIGURE 3. LOCUS OF ZEROS FOR INCREASING FULL 
ELASTIC COUPLING (R= I ), 8 = 1.0 ^ 
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FIGURE 4. 


ROOT LOCUS FOR INCREASING PITCH ANGLE. FULL 
ELASTIC COUPLING (R= I). VARIOUS cug. 
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FIGURE 6. LOCUS OF ZEROS FOR INCREASING PITCH ANGLE AS 
A FUNCTION OF R FOR =1.4. 
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* 1.4 

0 * 1.0 

FIGURE 7. ROOT LOCUS DEPARTURE ANGLES AS A FUNCTION 
OF R FOR SERIES AND PARALLEL MODELS 
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FIGURE 8. ROOT LOCUS FOR INCREASING PITCH ANGLE 
MATCHED STIFFNESS. COMPARISON OF SERIES 
AND PARALLEL MODELS. 
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FIGURE 9. COMPARISON OF ZERO LOCATION FOR RIGID MODEL 
AND FLEXIBLE BLADE MODEL AS A FUNCTION OF 
FLAP FREQUENCY AND INFLOW PARAMETERS. 
ELASTICALLY UNCOUPLED R=0. 
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FIGURE 10. LOCUS OF LEAD-LAG ROOTS . VARIABLE ELASTIC 
COUPLING. PARALLEL SPRING MODEL 
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FIGURE II. LOCUS OF LEAD- LAG ROOTS. VARIABLE 
ELASTIC COUPLING. PARALLEL SPRING 
MODEL 
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/ 

FIGURE 12. LOCUS OF LEAD- LAG ROOTS. VARIABLE ELASTIC 
COUPLING. PARALLEL SPRING MODEL 
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B BLADE PITCH ANGLE 


FIGURE 13. PRINCIPAL AXIS INCLINATION AS A FUNCTION 
OF BLADE PITCH ANGLE AND THE ELASTIC 
COUPLING PARAMETER R 
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APIEKDIX I 


FLAP-LAG EQUATIONS OF MOTION 

The nonlinear equations of motion for a variable property rotor 
blade undergoing flap^ lag and elastic twist deformations are developed 
in Reference h using a Rayleigh- Ritz procedure. Using a single mode 
for each type of deformation and Ignoring elastic twist altogether 
we may write the modal equations for flap-lag motion as follows. It 
has been assumed in developing the aerodynamic forces that the induced 
velocity is Independent of blade spanwise station. For convenience we 
have constructed nonlinear steady state (or static) equations and small 
perturbation (linear) dynamic equations about the steady state equilibrium. 
The latter are required for stability analysis. 

steady state 

o , vYv , „w9n , o „vw9 „y° 

v[M +K +K ]+wK = Qjf 

,pcw ^ o r^^Vw , .;ww9t , o .^vw9 „z° 

M +w[K +K ]+vK =Q 

PERTURBATION 

. nT.^UVW . O , . r^^UWWV O , 

v[M +K ]v+w[M w+M] 

+ v[M +K ,+ K ]+wK +M v^Q"^ 
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• rijW , 0 -, , ttz-Vw , 

v[M +K w]+w[K +K ] 


. ^ .,ww .. 

+ vK +M ■w=Q 


AERODYITAMIC FORCES 
STEADY STATE 


Q' 


o 


ac j-^ aX6v 2 „v , do ^2 .xxvt 

- p— tC! V. A - vf A + Q A } 

2 xn in a -* 


Q 


- pf {- Q= A^" + n a"3 


PERTURBATION 


Q' 


,y = 


ac r. rn «XbWV ^ „W-, , . (. „UW 0^ 

- p^ [w [Q A - 2v^^ A ] + V A + 2Q— A ]] 


Q 


ac r . _ .XWTT . . r r, 

p— Iw n A + V [Vj,^ A - 2QA ]} 


The several inertial, stiffness and aerod 3 mamic coefficient are given 
below, ijr^, are the lag and flap mode shapes assumed in the Rayleigh- 
Ritz analysis. Note that the mode shapes are dimensionless; however, 
()'=^(). 

MASS 


M 


.w 


R 


= - r m dr 

o ^ V 


M 


,w 


R 


= J m ilr^ dr 


M 


R 


= r m ilr® dr 
■Jo w 


M 


w 


R 


M =2QP rra'|f'Jf dr 
pc Jo ’ V ’^'W 
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P mr lit dr 
pc Jo w 

..uvw o n rR . , 

M = 2 m C V dr 

Jo ^v 

uwwv _ pR uww , 1 , , 

M =2Q inC V dr 

Jo V 

^Vw 

Jo '\r 

where 

, C“" = 20 jj m dn, C“* - -jj (♦;)= an, . -j^ (♦;)= an 

TVT 4. 4.1. 4. 74-VW ,,rU.VW ,,W ^ ,,WW ■ _2 

Note that K =-M •, M =-CrM , M =-0*^14 

STIFFMESS 

kV» =J«c“(t;)= dr 
c"“ = mT] dll 


w0 

K 

= [El / 

Jo z 

cos^ 0 

+ El / 

y 

sin® 

9] 

('P;')® dr 

,,ww9 

it 

= [El / 

Jo z 

'sin® 0 

+ El / 

y 

cos® 

9] 

(t ' " ) ® dr 

w ' 

k'™® 

- K*™® -J 

^ [El / 
o z 

- El / 

y 

-] sin 
■' 2 

20 

iV' V') dr 


0 is the rotation of the blade and in general consists of a constant part 

9 and a linear twist G, such that 
0 ^ 
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® - ®c R, 


In general, there -would be a segment of the blade/hub which is inboard 

of the blade pitch bearing and does not rotate with pitch. It is 

considered that this segment of the blade/hub system extends from the 

hub to a radius r^ and the stiffness characteristics of this segment 

are denoted by the subscript h. The outboard segment rotates with pitch 

and is denoted by the subscript b. In this case, therefore the stiffness 

coefficients are 
r. 


K EI^/ dr + (EI^/ cos^ e + El / sin^ 0 )(ilr'' 

h ^h a ^b, ^ 


dr 


K 


;ww0 


r 

= Jo ®V /r ® 

•^h ^-h ^ 


K 


,vw0 


= (El / - El ,) 
J r. z. y 


sin 20 


h ^b 


^") dr 

V ^w 


AERODYNAMIC 


^xev 


f 

‘t o 

r0 ij; 

^ V 

dr 



J O 

(It dr 

T Y 


aXXV 

= 

o 

r® \lr 
T-Y 

dr 

^XX0W 

= 

J 0 

r^ 0 dr 

^w 

A™ 

- 

J o 

r i!r dr 

X0WV 

= 

f 

J o 

r 0 i!/ 
^w 

■ij - 


,X0VW 
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.wv pR , , , .w 

A = dr = A 


„9w 

A 


.xwir 

A 



dr 


dr 

V 


If the same mode shapes are used for flap and lag deflections ('{f^ = ijr 
the following relationships hold among the coefficients: 



= = 

-Q^ 

^Vw 

^Vw ,,UVW 

= K = - M = - 

T,,UWWV 

M 

Vv 

K 

II 


Ti/r™ 

M 

po 



Define, noting that ( )' indicates a derivative with respect to x: 

M Ui) = ^ 

= !o^ 

M™ (xO = - 2 t dx 

d71 

K = sj' (V')^ -ax 

= Jo'' ^ 

X 


t) 
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If a segmented blade is considered with constant properties within each 
segment the several coefficients may be expressed more explicitly as given 
below. 

( 1 ) 

m = runnxng mass of segment located between 0 and x^. 

( 2 ) 

m = running mass of segment located between x^ and Xg. 


El 


( 2 ) _ 


stiffness of segment located between 0 and x^. 
stiffness of segment located between x^ and Xg. 


= pitch of segment located between 0 and x^. 

( 2 ) 

0p = pitch of segment located between x^ and Xg. 


= linear twist 


The tip of the blade is located at x^^. = 1. 

= m^^^ M (xi) + m^^^ [M (xg) - M (x^)] 

+ ... + m^^^ [M (xj^) - M 

R® I m^^^ (xi) 

+ m^2^ (xg) - (x,) 

+ • • • 

+ [sr -sr (x^_,)]} 
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^uvw ^ ^ Jl) ~uvW 


, (2) r«UVW I \ r,UVW / VT 

+ ^ [H (Xg) - M (xi)] 


+ ... 


, „(H) [i>^vvv 


STIFFNESS COEFFICIENTS 


K 


.w 6 2 


E® 

( 1 ) 


— |cos? K (xi) - [sin 2 0^^^] 01 I?" (xi)| 

( - 

- |si3 


El 

+ — ^— - |sin® 0^^^ K (Xi) + [sin 2 9^^^ 0i K (xj)] 
R 




El 


(M) , 

J,r" |cos= e(”) [i (Xg) - K 

- [sin 2 e, [?Tx^) - F (Xj,_^)]j 


ET<^) ( 

3in=9« 


[K (x^) - ic (x^_3^)] 


+ [sin 2 6i [? (x^) - i? (^. 3 ^)]| 


et 4) t 

K™® = —2— |sln= 9^^) K (xi) + [sin 2 0^^'] 0 


R® 

(1) 


El'" ( ri\ 

+ — 2L_ /cos® 0_ ^ K 
R® { 


X i? (Xx)| 
(xi) - [sin 2 0Q^b 01 (xi)| 


'f* 

' • « o 
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.vw® __w® 


K ■ = K 




2R" 


^ ^sin 2 K (xi) + 2 01 cos 2 9^^^ (xi)^ 


+ • • « 


=a=E 


+ . . . ' + 


2 


,m' ' — + m 


,( 1 ) 


m' ' X 


3 


F (Xi) 


(l) _ 

' rJ)xx / V 
“5 — ^ (xi) 


in 


(W) ~ ^-1^ 

2 


m 


(2)' 


(x,) - ie- (x,)] 


+ [ . . . 

(N) 

, 3 

2 

m"*) 


1 {...] 

t...) 


p (x^) - ? 

[?- (x^) - F- (X,.,)]. 


The various functions of mode shape are shown graphically in Figures A-1 
through A“7 for a Duncan Polynomial mode shape 

1 _ 2 ^ 1 A 

\^=2x - ^ x^ + X 

3 3 

The aerodynamic terms are as follows for this mode shape. 
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aerodiitamic coefficients 


11 

.Ur 

^xw 

A 

.289R3 

.XXV 

A = 

.225R® 

A^^ = 

.256R 

^xwv ^ 

•206R3 

0 w _ 

A “ 

• 206R ©i + .256 R 0^ 

X0WV 


A = 

• 172R^ 01+ .206 R® 0, 

.X0V 

A = 

.225R® 01 + .289R^ 6 

.XX0W 

A 

.I85R® 01 + .225R® 0 


The blade has a constant pitch 0^ and a linear twist represented by . 

To place these flexible blade equations in a form which can be conveniently 
compared to the rigid blade equations, the displacements are nondimensionalized 


by the rotor radius and each of the equations is multiplied by ( 


R 


I, Q 


0 


• 2 

where xs the flapping moment of inertia of the blade, = J ^ mr dr 


The time is nondimensionalized by the rotor RPM, n, such that 


V = 


R 


_o 
w ~ 


w 

R 


V = 


V 

QR ’ 


w = 


w 
Q R 


V = 


V It 
, W 


w 


R 


n® R 


The equations become, where the primes indicate derivatives of the mode 
shapes with respect to the dimensionless distance x, 
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7h 



w 


m 




-f m ^ ® dx 

O ^ V 


f m.x®dx 

n 


= - 111 


W 


ww 


m 


R® M™ 


r m ilf ^ dx 
Jo 


J m X® dx 


= - m 


ww 


w 


m 


wv 
= - ,m = 


R"M^ 


2 P r m il; * dx 
pc J o V w 


Ti Q 


L 


m X dx 


xw 


m 


RM^ 


P r m X ilr dx 
pc J o 


r m X® dx 
J o 


.uvw _ R^ 


m 


o P ^ , , 

2 m c <Ss dx 
Jo 


r m x^ dx 
J o 


where 


uw 


J? rtv')" 


. „a ,,u-wwv 

uwwv R M 


m 


I j Q 


2 ■[ ni dx 

Jo ^ t.t 


' w 


J m x^ dx 


Vw R-" K 


a 


I, Q 


r c^^(t < 3 x 

r m x^ dx 

J o 


mv ^ ^ , __ 

c - 2 m dT1 
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and 


, Vw 
k 




Jo 


mv 


(f ) 


8 


J 


m X 


dx 


dx. 


, Vw uvw 
k = -m 

The aerodynamic terms are, denoting by small a*s the dimensionless form 
of the A’s times 4 , 


X0V 4 A' 


xQv 




= 4 X 0 dx 


and X = 


Q R 


8Q'' 


pac R® Q® 


r , X0V - 3 V , do „XXVo 

-[-Xa -Xa+ a J 

a 


8q 


pac R® 


r XX0W , XWn 

= - [-a - X a J 


8Q' 


,y 


pac R^ Q® 


r-L ( X0WV , „ 1 w^ , ^ ^ 9w , ^ do xwvv-, 

-[w(a +2Xa ) + v (-X a +2 a )] 

a 


8q" 


pac R® 


r— XWV , J. / , W ^ X0VWs, 

[wa +v(-Xa -2a 


The stiffness terms are 


= k'''' Jo 


Vv R 


J m X® dx 


•ir 

^Vw _ R K ^ J_0 

fi® J m X® dx 
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vt/6 K 


3 






f (EI / cos^ 0 + El /sin^ 0)(’!f^^)^ dx 
0 o z y 


y 


f m X® dx 
<J Q 


k 


WW0 R*^ K 


3 T^WwS 


Ii Q' 


R® 


(EI^/ sin® 0 + Ely/ cos® dx 

r in X® dx 
J o 


k 


VW0 R"' K 


2 ..we 


Q® 


f" (EI EI /) — - g-— ijf") dx 

o ^ z Y 2 ' 


R® Q® 


r m X® dx 
J o 


The equations become, using this notation, 

Steady State Equations 

C 

1 V „o , , VW0 -o y X0V . . 2 V do „xxv. 

(m +k +k )v+k ’w=*(Aa +X a- a } 

xw , Vw , , ww9 V _o , , w6 „o 7 / xxBw , . xw. 
m +(k +k )w+k v = '^(a +Aa; 


Perturbation Equations 

— / uww _o , W-Vx , _ / w , , V*v , , W0 \ 
■w(m w+m)+v(m +k +k ) 


, _ , VW0 ,X W Y / . / X0WV ^ WYx . — , 0W do „XWVsv 

+ wk +vm = ^ (“W (a + 2 a a ) + v (+A a - 2 a J) 


— / w . , Vw _o X . _ Vw , , wwe X _ . vwQ X - 
v(m +k •w)-s-w(k +k )+v(k ) 

“ w Y / — XWV . w , o X0Wxx 

+ wm =* (-wa +v(+Xa +2a )) 
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Dividing the perturbation equations by the coefficients of the 


acceleration terms and collecting terms on the left hand side. 


-• V / ^do ,a^ 


0 w \ . / vv , , Vv , - w0 

/-a \ ./111 +k 

(— jv- ) - A. (-U ) 1 V + — | V 

w ' ' w / V vv 

m m / ' m 


. - im 
+ w 


WJ0 


UWW W 3C0WV WV \ 

w- - + - 5 (V - ™ (V) = 0 

m m m ra ' m 


^2CWV ^ ww0\ / vw , Vvw 

-,Y/a NJ-,(k +k ,/k ^_o 

o w \ w / I' ww ' w ' 

m ' m / \m m 


I 


vw X0VW \ , VW0 

(1^ + 2 ) + ? (K^) = 0 

' ww ww V ' ww 

m m / m 


For zero precone. 


m =0 


vw ^ 
m =0 


and defining the dimensionless uncoupled frequencies as 
Vw , , WW0 


2 k" ' + k 
P = 


ww 


m 


w , , Vv . , w0 
„2 m + k + k 

q = j; 

w 

m 


the coefficients of the Coriolis terms as, 


-w m 


uwwv 


w 


m 


..ww k 


Vw 


ww 


m 
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and the elastic coupling terms as 




vw9 


w 


m 




vw6 


ww 


m 


the equations are 


a 




C xwv 0 w 
do /a \ , a 

) — A — ^ 

^ tj- 


W 


m 


w 


m 


V + q^ V 


+ 


. / ,, xQ^wv wv \ 

— _o , Y /a , a \\ , „ 

W I C W + (-TS — + 2X -^r— ) + C-, w = 0 

V o o w w^l -*■ 

' m m / 


. xmr / a, "vw X0W \ 

w + ^ (- 1 ; — ) w + p w + V (C w - -Q (,A-ii— -i- 2 -77“ — M 
o ww ' \ o o ww ww / 

m mm 


\ 1 „WVf — ^ 

)/+ Cl V = 0 


If the mode shapes are assumed to be the same for both deflections then, 


ww w 
-m = m 


✓'WV ww 

c = c 

000 


r. _ _ pWW 

“^1 ~ '-'S 


and defining 


a 

a = 


IX 


11 


w ww 
m m 


the equations are, recalling that X is the inflow angle at the blade tip, 
C 

^ , Y /o do ~xwv i ~0w,, J. , 3 — 

V + -k (2 a -Aa )v+qv 

o a 


+ ^ (-C^ ^ 2 X a^) ) + Ci w =0 
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w 


+ y- 


~XViTV — 

a w 


p® w 


+ V [C 


_o 

■W 




a 


+ 2 


~x6wO\ 
a )j 


+ Cl V = 0 


Now these equations are in a form which may be compared with the rigid blade 
equations. If a rigid blade model is employed, then the spring constants for 
the rigid blade model are identified in such a way that the parallel spring 
model gives identical coefficients to the stiffness terms represented by 
p^, q® and Cj as well as in the calculation of the coning angle w^. The 
significance of these changes is considered elsewhere in this report. 

The steady state equations are 

2—0 ^ — o y (\ ~x0v ,2 

q V + CiW = g lla +X^a 

2 — o , „ — o Y /~xx0w , -> ~xw\ 

p w +CjV = g (a +Aa I 

The rigid blade terms can be obtained by substitution of ^If = x in 
the expressions given above. 
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TABLE I 


COMPARISON OF 

COEFFICIENTS IN FLAP-LAG EQUATIONS FOR RIGID BLADE 
AND DUNCAN POLYNOMIAL MODE SHAPE. UNIFORM MASS DISTRIBUTION 

Aerodynamic Terms 



4 = ill = X 


x 9 v 

a 

1.333 

1.156 0^ + .900 0 1 

V 

a 

2.00 

1.600 

xxv 

a 

1.00 

.900 

XX0W 

a 

1.00 0 + .800 0., 
o ^ 

.900 0^ + .738 61 

xw 

a 

1.333 

1.156 

X 0 WV 

a 

1.00 0 + .800 0n 

o ^ . 

.824 0 + .687 01 

0 _ 

w wv 

a = a 

1.333 

1.028 

0W 

a 

1.333 + 1.00 0 i 

1.(^8 0 + .824 0 1 
0 ^ 

xwv 

a 

1.000 

<M 

CO 

Mass Terms (uniform mass blade) 





i(r = ,1| = X 

V w 

^ ^ = 2x^ - 1 x=^ + 1 X - 

V V 3 3 

w 

m 

- 1 

- .770 

w 

m 

1 

.770 

w 

m 

pc 

1.542 

xw 

m 

pc 

•867 Bp, 

UWWTT 

-m 

2 . 

1.625 


8l 


COlH U>|H 












Stiffness (uniform stiffness blade) 



^ = X 

t = 'll = 2x® - ~ 1 

W 3 3 

^VVT 

2 

1.625 


2 

1.625 

11 

1 

. 90 ^. 
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X = r /R 

FIGURE A-1 . MODE SHAPE INTEGRAL 
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4,0 



FIGURE A-4. MODE SHAPE INTEGRAL 
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